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ABSTRACT 

A previously unseen maser in the J = 1 - transition of HCO + has been detected by the Combined Array for 
Millimeter-wave Astronomy (CARMA). A sub-arcsecond map was produced of the 2 arcmin 2 region around 
DR21(OH), which has had previous detections of OH and methanol masers. This new object has remained un- 
detected until now due to its extremely compact size. The object has a brightness temperature of > 2500 
K and a FWHM linewidth of 0.497 km s" 1 , both of which suggest non-thermal line emission consistent 
with an unsatur ated maser. This object coin cides in position and velocity with the methanol maser named 
DR21(OH)-l by lPlambeck & Mentenl (0,990). No compact HCO + emission was present in the CARMA data 
to wards the other methanol m asers described in that paper. These new results support the theory introduced 
in Plambeck & Menten ( 1990) that these masers likely arise from strong outflows interacting with low mass, 
high density pocke ts of molecular gas. This is further supported by recent observations of a CO outflow by 
Zapata et al. (2012) that traces the outflow edges and confirms that the maser position lies along the edge of the 
outflow where interaction with molecular tracers can occur. 
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1. INTRODUCTION 

DR21(OH) is a part of a massive star forming region in 
the DR21 molecular cloud complex. It lies along a ridge 
of strong molecular emission extending from DR21 in the 
south to W75N in the north. As with many objects of its 
type, DR21(OH) also contains numerous sites of strong maser 
emission in many species. Named for its particularly strong 
OH-masers, it was determined that the OH and H2O masers 
were strongly co incident with millimeter continuum sources 
dPadin et alj|1989h . 

DR21(OH) is particularly significant due to the detec- 
tion of methanol maser s over a wide range of frequencies. 
IB atria & Mentenl (fl98 8) detected metha nol masers in the 8 1 .6 
GHz a nd 84 GHz methanol transitions. Plambec k & Mentenl 
( 1990) conducted a study of 95 GHz methanol emission and 
found four strong methanol masers across a region associ- 
ated with the MM2 millimeter source. ISlvsh et alTd 1 9971) pre- 
sented detections of a previously u n seen 133 GHz methanol 
maser. More recently. lAraya et aD (120091) present an exten- 
sive survey of over 30 44 GHz methanol masers, which ap- 
pear to trace the shock fronts of two bow shocks along the 
red and blue shi fted lobes of an outflow generated by MM2. 
iFish et ail d201 ll) present an analysis of all methanol maser 
detections and determine that all of the known masers appear 
to come from the interface region along the shock fronts of the 
outflow. The brightest of these masers appear along the west- 
ern tip of the outflow and occur in a narrow velocity range of 
approximately 0.3-0.5 km s -1 . Additionally, they determine 
that all the methanol maser transitions are Class I masers, 
which arise from collisional excitations and are considered 
to be caused by shocks, especially those that arise from out- 
flows, as opposed to Class II masers, in which the pumping 
mechanism is primarily derived from external radiation. Re- 
cent submillimeter c ontinuum and 1 mm CO detections of the 
outflow morphology (Za pata et alj |2012) give further support 
to the above model. In this letter we report the first detection 
of an HCO + maser in the interstellar medium. 

2. OBSERVATIONS 

DR21(OH ) was observed as part of an ongoing survey of 
sources (see lHakobian & Crutcherl 1201 ll) with strong mag- 



netic field detections dFalgarone etalj |2008; Crut cherl [19991) . 
A 2 square arc minute region around DR21 (OH) was mapped 
with observations by CARMA dBock et al]l2006l) in the C, D, 
and E array configurations with a combined effective beam 
size of 6.3" x 4.0", and a spectral resolution of 0.41 km 
s -1 . This was achieved by observing in the 8 MHz band- 
width correlator mode with 63 channels. The J = 1-0 tran- 
sition of HCO + , with a rest frequency of 89.188526 GHz, 
was simultaneously observed with HCN, N2H + , and two con- 
tinuum bands. A bright, unresolved, and unusually narrow 
line feature was detected in the map of HCO + (Fig. [U, at 
20 h 38 m 59.3 s +42°22'49.0" (J2000) with a peak at V LS r = 
0.78 km s" 1 . Each track was inspected to rule out the pos- 
sibility of a transient instrumental issue; however, the com- 
pact source was visible in each track. The spatial and spectral 
location of the source additionally does not change between 
each of the tracks, indicating that an instrumental issue is un- 
likely. Follow-up observations in CARMA B-array configu- 
ration (Fig. O were performed in December 2011 in order 
to further constrain the physical size and brightness of the 
source. In total, five tracks of data were observed for a to- 
tal of 18.75 hours: 4 hours each in B, C, and E arrays, and 
6.75 hours in D array (over two tracks). Data were calibrated 
and imaged with the CARMA MIRIAD software package. 
The phase calibrator used in all tracks was 2038+513. Due 
to the strength of the calibrator, we were able to use the stan- 
dard phase calibration techniques. The B-array track has an 
improved spectral resolution of 0.14 km s" 1 over the other 
tracks, due to being observed with the expanded CARMA 
8-band correlator in the 8 MHz bandwidth mode with 384 
Hanning smoothed spectral channels (Figure|2]shows the 192 
spectrally independent channels). 

3. data / analysis 

Using the high resolution B-array data, we find that the an- 
gular size of this object is less than ~ 0.8 arcseconds and that 
it has a brightness temperature > 1900 K. Since the source 
remained unresolved in B-array, the parameters of this object 
can be further constrained by using only long baseline com- 
ponents of the data set. This effectively "resolves out" larger 
scale structure, leaving a constrained map behind. To achieve 
this, we limited the map to use baselines larger than 140 kA 
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Figure 1. Comparison of maser emission in CARMA C, D, and E-arrays. Spectra (right) are of the peak positions in each map (left) of the 0.78 km s 1 channel, 
visible in shaded contours. The y-axis of the spectra are shown in both units of Flux density (Jy beam -1 ) and Brightness Temperature (K) to emphasize the effect 
of the decreasing source/beam size. The beam sizes are 2.0" x 1.5", 6.4" x 3.7", and 1 1.9" x 8.4" for the C, D, and E-array maps, respectively. The contour levels 
for the C-array map are 0.25, 0.50, and 0.75 times the peak flux density. The contour levels for the D and E-array maps are 0.01, 0.075, 0.15, 0.30, 0.45, and 0.7 
times the peak flux density. The RMS noise of the spectra are 0.042, 0.028, and 0.056 Jy beam -1 for the C, D, and E array maps, respectively. 



(~ 470 meters). This procedure resulted in a map with an 
effective beam size of ~0.54" x 0.52" and a brightness tem- 
perature of > 2500 K. The source appears to be unresolved 
even at this resolution (Fig [2] top panel). 

If the source were dominated by thermal emission, a lower 
bound on its linewidth can be estimated by calculating the 
degree by which the emission would be thermally broad- 
ened. The calculated full-width at half maximum (FWHM) 
linewidth would then be purely a function of the effective 
brightness temperature: 

/8Jtrin2 A/ Av 
A/ = \/ =— /o, — = — , m = Nm p 

V mc A jo c 

= 0.2139 J— kms" 1 
V N 

where T is the brightness temperature of the source and N 
is the atomic weight of the molecule. The factor of 2\/21n2 
comes from the conversion of a Gaussian standard deviation 
to FWHM. For a source with a brightness temperature of 2546 
K and RMS noise of 8.6 K (0.015 Jy beam -1 ), the thermal 
linewidth would be 2.077 ± 0.004 km s" 1 , which is signifi- 
cantly greater than the observed linewidth of 0.497 ± 0.002 



km s (obtained by least-squares gaussian fitting to the line 
profile). The maser linewidth corrected for the instrumen- 
tal smoothing is 0.477 ± 0.002 km s" 1 . The observed high 
brightness temperature and linewidth narrowing are good in- 
dicators that this object is dominated by non-thermal emis- 
sion. An assumption as to the type of maser emission can 
be made from this linewidth analysis. Saturated masers can 
have linewidths up to the thermal linewidth, while unsatu- 
rated masers will have a linewi dth narrower than the thermal 
linewidth by a factor of 4 to 5 dReid & Moranll 19881) . Since 
our object has a linewidth that is 4.2 times narrower than the 
th ermal linewidth, it is c onsist ent with an unsaturated maser. 

Plambeck & Menten (1990) observed this region with 
BIMA at the 95 GHz methanol line and discovered four 
methanol masers connected by large scale methanol emission. 
The brightest of these four sources, DR21(OH)-l, is centered 
at 20 h 38 m 59.24 s +42°22'49.04" (J2000) with a peak at V LS r= 
0.32 km s" 1 , which is approximately 0.1" from the center of 
our measured HCO + peak. The HCO + position has an ap- 
proximate positional error of 0.2", and the methanol positions 
have an error of 0.3". From these data they estimated their 
methanol maser had an angular diameter <4.4" and a bright- 
ness temperature > 760 K. Three other methanol masers were 
also observed; however, there is no evidence of companion 
HCO + masers (Fig. [3}. This could possibly be due to their 
relatively weaker strength or due to the fact that the condi- 
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Figure 2. Comparison of maser emission in CARMA B-array. Spectra (right) are of the peak positions in each map (left) of the 0.78 km s~' channel, visible in 
shaded contours. The y-axis of the spectra are shown in both units of Flux density (Jy beam -1 ) and Brightness Temperature (K) to emphasize the effect of the 
decreasing source/beam size. The maps feature an inlay showing the 4 arcsecond 2 region around the maser due to its compact size. The top panel is of B-array 
data with baselines > 140kA. The beam sizes are 0.54" x 0.52" for the top panel, and 0.85" x 0.70" for the bottom panel. The contour levels are 0.25, 0.50, and 
0.75 times the peak flux density. The RMS noise of the spectra are 0.015 and 0.026 Jy beam -1 , for the top and bottom map, respectively. 



tions necessary to produce an HCO + maser do not exist at 
these other positions. The co-location of the HCO + object and 
the met hanol maser suggest that bo th arise from similar con- 
ditions. Plambeck& Menten (1990) hypothesized that such 
regions could be created from the interaction of an outflow 
and small clumps of dense molecular gas. 

In comparing the B-array maps (Figure 01, it appears as if 
the maser could be partially resolved since its intensity drops 
from 7.4 Jy beam" 1 to 4.7 Jy beam" 1 between them, with 
slight N-S structure. Ara va et ail (|2009) reports that there are 
three methanol masers within 2" of the HCO + emission; one 
strong and two weaker masers, slightly north and south of the 
strongest maser. If the HCO + maser is also a complex of three 
individual sources oriented in a N-S direction, the highest res- 
olution B-array map may be beginning to distinguish them, 
giving an appearance of being resolved. In the D and E ar- 
rays, the sharp maser "peak" extends ~ 7.5 Jy beam" 1 above 
the extended emission, in agreement with the flux density in 
the B-array spectrum. 

A small amount of HCO + absorption is visible in the E- 
array map; however, it is not apparent that the maser linewidth 
is affected by self-absorbtion. In this map, the shortest base- 
lines will sample large enough size scales to include the west- 
ern extension of the continuum that peaks with DR21(OH)- 
MM2. Furthermore, the HCO + emission includes two veloc- 
ity components associated with extended gas, one at -4.7 km 
s , and another peaking at ~ km s" 1 . The peak positions of 
these two velocity components correspond to the continuum 
peaks of MM1 and MM2 respectively. The high velocity wing 
of the -4.7 km s" 1 component is affected by absorption by the 
continuum in the E-array map, however, the km s" 1 compo- 
nent is not. The D and C array maps with smaller beam sizes 
(longer baselines) are not affected by this continuum contam- 
in ation or the extende d HCO + emission which is resolved out. 

Zap ata et al.l d2012) performed 1mm observations of several 



spectral lines around DR21(OH). Included with these obser- 
vations are CO(2-l) observations which trace an outflow from 
DR21(OH)-MM2. This outflow appears to be in the plane of 
the sky (Fig. [3). From this figure, we can see that both the 
methanol maser and our HCO + object appear along the edge 
of the outflow. This region would be highly shocked, and 
the energy released from this interaction has the potential to 
p ower the maser. 

iGoldsmithl (119721) analyzed the J = 1 - transition of CO 
and determined that a large range of rotational excitation tem- 
peratures, including population inversion, can be produced 
through collisional excitation in the range of kinetic temper- 
atures and densities found in molecular clouds. It was also 
concluded that other linear molecules with simple rotational 
structure, such as HCO + , would have the same result. Us- 
ing th e RADEX radiative transfer package (IVan der Tak et al.l 
120071) . we performed a test calculation that showed that for 
collisional interaction such as is suggested here in DR21(OH), 
it is possible to achieve population inversion in the HCO + J 
= 1-0 transition with gas densities of ~ 5 x 10 s cm" 3 and 
kinetic temperatures ^50 K (consistent with kinetic tempera- 
tures in outflow shocked regions). This result further supports 
the conclusion that this HCO + is a maser. 

4. CONCLUSION 

Observations of DR21(OH) have revealed the presence of a 
compact object which is dominated by non-thermal emission. 
The extremely compact size of this object coupled with its 
large brightness leads to the following conclusions: 

• The source is co-located with a known strong methanol 
maser. 

• It lies along the edge of an outflow which gives sup- 
port to previous theories that masers can arise due to the 
interactions of high velocity outflows with cold, dense 
clumps of molecular gas. 

• Its small spatial size most likely prevented its detection 
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Figure 3. Diagram of the position of the HCO + maser with respect to o ther features of DR 21(OH). The shaded contours are of the 112 GHz continuum of 
DR21(OH). The dotted arrows represent the CO outflow as observed by Zapata et al. 12012). The shaded oval represents the extent of the HCO + object as 
reported in this paper. The stars represent the positions of methanol masers from Plambeck & Menten 1 1990). The black curves represent the two sets of bow 
shocks determined by the loci of 44 GHz methanol masers from Araya et al. 1 20Q3). The positions of the masers along the edge of the outflow give strong support 
to the theory that shocked interactions between outflows and small, high density clumps of molecular gas give rise to these objects. The contour levels of the 
continuum emission are 0.3, 0.4, 0.49, 0.54, 0.7, 0.9, and 0.99 times the peak flux density of 0.103 Jy beam -1 in order to show the peaks of both continuum 
sources. 



before now; emission from this source would be beam 
diluted to levels indistinguishable from thermal emis- 
sion. 

• This object is very likely an unsaturated maser, the first 
observed in HCO + . 
In order to further confirm that this source is indeed a maser, 
future observations to look for anomalous level populations in 
higher order HCO + transitions should be performed. 
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